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ANALYSIS OF BROKEN PLANE TRAJECTORIES FOR THE PIONEER -F MISSION 

by Fred Teren 
Lewis Research Center 

SUMMARY 

A study was conducted to determine if broken plane trajectories can be used to in- 
crease the duration of the Pioneer- F launch opportunity, or, equivalently, to increase 
payload capability for a fixed launch opportunity duration. Optimum single conic and 
broken plane trajectories are calculated for the launch day - trip time combinations, 
and minimum trip time trajectories (for a given payload mass) are determined for each 
launch date. The maximum mass of spacecraft propellant required for the broken plane 
trajectories is calculated. Then, new minimum trip time trajectories are calculated in 
which the added spacecraft propellant mass is optimally utilized for all launch days. 
Launch opportunity duration is presented as a function of spacecraft mass for the single 
conic and broken plane opportunities. Values of launch vehicle injection hyperbolic 
velocity and declination are presented, along with the values of optimum spacecraft pro- 
pellant mass and correction time. The results show that the use of broken plane trajec- 
tories can provide a small increase in launch-opportunity duration and/or spacecraft 
mass. 


INTRODUCTION 

The Pioneer- F is a Jupiter flyby mission to be launched in 1972 on an Atlas-Centaur 
launch vehicle with a solid- propellant third stage (TE 364-4). As has been typical for 
interplanetary missions, preliminary trajectory planning called for single conic trajec- 
tories to be used exclusively in establishing the launch opportunity; that is, all the re- 
quired spacecraft velocity is added by the launch vehicle in the vicinity of the Earth. 

The heliocentric trans-Jupiter trajectory then consists of a Keplerian arc (ellipse) be- 
tween Earth and Jupiter, modified only by midcourse corrections performed by the 
spacecraft to correct for velocity errors resulting from launch vehicle guidance system 


errors. 



The present study was conducted to determine if broken plane trajectories can be 
used to lengthen the Pioneer- F launch opportunity and/or to increase spacecraft mass 
capability. With broken plane trajectories, the spacecraft midcourse-correction pro- 
pulsion system is used to change the trajectory velocity vector in heliocentric space. 

For certain launch day - arrival day combinations, usually when the heliocentric travel 
angle is near 180 degrees, such a maneuver can increase payload mass capability by 
reducing the launch vehicle velocity requirements. 

Broken plane trajectories have been shown by several authors to be useful in re- 
ducing energy requirements for interplanetary missions. A recent paper by Gobetz and 
Doll (ref. 1) presents results for Mars orbiter missions, with primary emphasis on the 
1971 and 1973 opportunities. However, reference 1 and earlier papers (refs. 2 and 3) 
used velocity increment, AV, as the performance index. In actual mission applications, 
it is desired to maximize either payload capability or launch opportunity duration, sub- 
ject to mission constraints. This requires consideration of spacecraft mass and pro- 
pulsion system characteristics, and of launch vehicle injected mass capability as a 
function of magnitude and direction of the geocentric hyperbolic velocity vector. 

Results are presented for spacecraft mass capability as a function of trip time for 
the various launch dates for single conic and broken plane trajectories. These data are 
crossplotted to obtain minimum trip time as a function of launch date for selected values 
of spacecraft mass. The number of possible launch days is then determined as a function 
of spacecraft mass for different values of maximum allowable trip time. In addition, 
the values of hyperbolic velocity, declination of the geocentric asymptote, optimum 
spacecraft propellant mass, and time of spacecraft propulsive maneuver are also pre- 
sented. 


METHOD OF ANALYSIS 

The mission strategy assumed in this study is to select, for each launch date, the 
shortest trip time consistent with launch vehicle capability and spacecraft mass. Launch 
vehicle injected mass capability is calculated as a function of hyperbolic excess velocity 
and the declination of the geocentric asymptote. The spacecraft propellant mass re- 
quired for the plane change is calculated by using the specific impulse of the spacecraft 
midcourse correction propulsion system, which was initially designed to correct for 
guidance system errors. To use this spacecraft propulsion system to effect plane 
changes, additional propellants, propellant tank capacity, and possible structural modi- 
fications would be required for the spacecraft. The existing spacecraft propulsion 

system uses monopropellant hydrazine, with a specific impulse I of 225 seconds. 

sp 
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This value of I is used in this report. However, the increase in propellant tank 
S P 

mass and, possibly in structural mass, is not evaluated in this report. 


Calculation of Launch Opportunity 

The launch opportunity is determined by calculating, for each launch day - trip time 
combination, the maximum spacecraft mass which can be delivered, both for single conic 
(no spacecraft propulsive maneuver used) and broken plane trajectories. Then, for 
each launch day, the minimum trip time is determined (for a given spacecraft mass). 
This trip may be either a single conic or a broken plane transfer. The required space- 
craft propellant loading is the maximum required on any day in the launch opportunity. 

If it is assumed that this spacecraft propellant must be tanked for all days in the oppor- 
tunity (including those days when single conic trajectories are optimum), a new set of 
optimum transfers must be determined, in which the full amount of added spacecraft 
propellant must be consumed in an optimum manner on each launch day. The resulting 
minimum trip times for the modified opportunity are longer than the previously calcu- 
lated values, since it was not optimum to use the full amount of propellant except on one 
launch day. 

The Fletcher- Powell gradient method (ref. 4) is used to determine the mass of 
spacecraft propellant, and the time and direction of spacecraft burn, which maximize 
payload capability for each launch day - trip time combination. 


Calculation of Trajectories 

The single conic launch opportunity is determined by calculating the spacecraft mass 
capability for each launch day - trip time combination. The required trajectories are 
calculated by solving Lambert's problem (ref. 5). Then, the required AV at Earth is 
rotated from heliocentric to geocentric coordinates, and the required hyperbolic velocity 
V H and declination 6 and right ascension <p of the outgoing asymptote are determined. 
Launch vehicle injected mass capability is calculated as a function of Vjj and 6 by 
using an empirical formula obtained from integrated launch vehicle trajectories. The 
value of (p simply determines the optimum launch time, corresponding to a zero- 
duration launch window. However, since a finite launch window is required on each 
launch day, the injected weight must be reduced below its value for the optimum launch 
time and cp. The reduction is not considered herein, but unpublished results have 
shown the mass loss to be from 5 to 10 pounds mass (2. 3 to 4. 6 kg) for a 30-minute 
launch window . 
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The optimality of the single conic trajectories is checked by calculating the magni- 
tude of the variational primer vector along the reference trajectory. This procedure is 
described in appendix B. (All symbols are defined in appendix A. ) If the single-conic 
trajectory is found to be nonoptimum, an optimum broken plane trajectory is calculated 
for the same launch day - trip time combination. The optimum broken plane trajectory 
is determined by using the Fletcher-Powell gradient method (ref. 4) to optimize the 
time, magnitude, and direction of the spacecraft impulse. 

Once the optimum trajectories are determined, the launch opportunity is defined by 
selecting, for each launch day, the shortest trip time which allows the required space- 
craft mass to be achieved. The resulting trips may be either single conic or broken 
plane transfers. 


RESULTS AND DISCUSSION 


The spacecraft mass capability is shown as a function of trip time for three launch 
days in figure 1. Net spacecraft mass is defined as the injected spacecraft mass less 



Figure 1. - Maximum net spacecraft mass as a function of trip time for three launch days. Single 
conic and broken plane trajectories. 
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the amount of spacecraft propellant added for the spacecraft impulse (for the broken 
plane trajectories). For launch on February 23, 1972, which is early in the launch 
opportunity (as will be shown later), broken plane trajectories are optimum for all trip 
times and all values of spacecraft mass. However, the difference between single conic 
and broken plane capability is small for trip times up to 650 days. For a spacecraft 
mass of 500 pounds mass (227 kg), single conic trajectories could be used with a very 
slight increase in trip time. However, for a spacecraft mass of 510 pounds mass 
(232 kg) or greater, single conic trajectories cannot be used. Broken plane trajectories, 
on the other hand, can be used for values of spacecraft mass up to 535 pounds mass 
(243 kg) with trip times of 850 days or less. 

For launch on March 6, 1972 (the center of the opportunity), single-conic trajec- 
tories are optimum for all values of spacecraft mass, assuming that the shortest trip 
time is desired. The single conic trajectories are also optimum for a March 18, 1972, 
launch date (near the end of the single conic opportunity) if the spacecraft mass is 
520 pounds mass (236 kg) or less. If the spacecraft mass is between 520 and 575 pounds 
mass (236 and 261 kg), broken plane trajectories may be used. The required trip times, 
however, are rather long, ranging from 850 days for a 520- pound-mass (236-kg) space- 
craft to more than 1100 days for a 570- pound-mass (259-kg) spacecraft. 

Similar data are displayed in figure 2 for a family of launch dates. Figure 2 differs 



460' 1 I l I I 1 I I I I I I I I 

500 600 700 800 900 1000 1100 1200 


Trip time, days 

Figure 2. - Maximum net spacecraft mass as a function of trip time for a family of launch dates. 
Single conic and broken plane trajectories (only optimum shown). 
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from figure 1 in that only the optimum trajectory (single conic or broken plane) is shown 
in figure 2 for each launch day and trip time. The following trends are readily observ- 
able from figure 2: For early launch dates, broken plane trajectories are optimum. As 
the launch date increases to about February 27, 1972, single conic trajectories become 
optimum and remain so up to about March 18. At this time, broken plane trajectories 
become optimum again. 

The data in figure 2, and other data of the same type for other' launch days, are 
crossplotted to determine the launch opportunity. For each launch day, the minimum 
trip time is determined for a given spacecraft mass, and plotted as a function of launch 
date in figure 3. The required trip times increase with spacecraft mass. For space- 
craft masses of 540 pounds mass (245 kg) or less, broken plane trajectories supply 
additional launch days prior to the opening and after the close of the single conic oppor- 
tunity. The required broken plane trip times for launch in late February 1972, are in 
the same range as those required for single conic trajectories. However, for broken 
plane trajectories in late March, the trip times increase rapidly. The curves were 
arbitrarily terminated at a trip time of 1200 days, although they could have been con- 



Launch date (1972) 


Figure 3. - Minimum trip time as a function of launch date. 
Optimum spacecraft propellant mass. 
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tinued to give additional launch dates with trip times longer than 1200 days. For space- 
craft mass greater than 540 pounds mass (245 kg), broken plane trajectories do not exist 
prior to the opening of the single conic opportunity. 

From figure 3, the number of available launch days may be calculated for each 
value of spacecraft mass. For a spacecraft mass of 520 pounds mass (227 kg), a 23-day 
launch opportunity is available with single conic trajectories. Maximum trip times of 
725 to 875 days are required. The launch opportunity can be increased to 27 days by 
using broken plane trajectories with no increase in maximum trip time, or to 39 days 
if 1200 day trips are permitted. If it is decided that a 23-day launch opportunity is 
adequate, then broken plane trajectories allow an increase in spacecraft mass to 
540 pounds mass (245 kg) for a 975-day maximum trip time, or to 560 pounds mass 
(254 kg) if 1125-day trips are allowed. Thus, it is seen that the broken plane trajec- 
tories can be used to increase the number of available launch days, or the spacecraft 
mass for a fixed number of launch days, or a combination of both. 

The mass of spacecraft propellant required for the broken plane transfers is pre- 
sented in figure 4 as a function of launch date, for various values of spacecraft mass. 



Q l I I I I I I I I I II I I I I I 
2/19 2/23 2/27 3/2 3/6 3/10 3/14 3/18 3/22 3/26 3/30 4/3 

Launch date (1972) 


Figure 4. - Optimum spacecraft propellant mass as a function of launch date for 
broken plane trajectories. 


7 





2/19 2/23 2/27 3/2 3/6 3/10 3/14 3/18 3/22 3/26 3/30 4/3 


Launch date (1972) 

Figure 5. - Optimum time of spacecraft correction as a function 
of launch date for broken plane trajectories. 

The corresponding time of the spacecraft propulsive maneuver is presented in figure 5. 
These values correspond to the minimum trip times for each launch date and spacecraft 
mass, as presented in figure 3. It can be seen from figure 4 that about 100 pounds mass 
(45 kg) of spacecraft propellant is required to take full advantage of the broken plane 
trajectories. Increased propellant tank capacity (hence, mass) would be required for 
this additional propellant mass. In addition, some structural modifications might be 
required for the spacecraft. These mass increases are not considered in this report. 
Figure 5 shows that the optimum spacecraft burn times occur about 100 to 200 days 
after launch for broken plane transfers early in the opportunity, while the corresponding 
times late in the opportunity range from 250 to 350 days after launch. 

The required values of launch vehicle hyperbolic velocity and declination of the geo- 
centric asymptote are presented in figures 6 and 7, respectively. For single conic 
transfers, the hyperbolic velocity is very nearly constant for a given value of spacecraft 
mass. The required launch-vehicle hyperbolic velocity decreases when broken plane 
transfers are used, because the spacecraft supplies some of the required velocity. The 
values of declination shown in figure 7 become less negative when broken plane transfers 
are used. This is reflected in less southerly launch azimuths. 

From the results presented in figure 4, about 100 pounds mass (45 kg) of spacecraft 
propellant mass should be added to take full advantage of the broken plane trajectories. 

A question that arises is how this propellant should be utilized in the single conic oppor- 
tunity and on those days when broken plane trajectories do not require the full 100 pounds 
mass (45 kg). The optimum procedure would be simply to load the optimum amount of 
spacecraft propellant (if any) each launch day. If this procedure is followed, the launch 
opportunity will be unaffected, except for a slight increase in spacecraft system mass 
resulting from increased tankage and structural modifications. However, it is possible 
that spacecraft propellants may not be varied throughout the launch opportunity. In this 
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Figure 6. - Injection hyperbolic velocity as a function of launch date. 
Optimum spacecraft propellant mass; single conic and broken plane 
trajectories. 



Launch date (1972) 

Figure 7. - Declination of the geometric asymptote as a func- 
tion of launch date. Optimum spacecraft propellant mass; 
single conic and broken plane trajectories. 







case, a fixed spacecraft propellant load must be tanked and utilized in an optimum 
manner. 

If the spacecraft is forced to contain propellants in excess of the optimum amount, 
the spacecraft mass capability will decrease for each launch day - trip time combination. 
The trip times shown in figure 3 cannot then be used. A new set of optimum trajectories 
must be obtained for the launch day - trip time combinations, in which 100 pounds mass 
(45 kg) of spacecraft propellant must be optimally utilized for each case. 

Spacecraft mass capability is presented as a function of trip time for several launch 
days in figure 8 for the restricted optimum trajectories. Comparison with figures 1 
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Figure 8. - Maximum net spacecraft mass as a function of trip time for several launch 
days. Spacecraft propellant mass, 100 pounds mass (45 kg). 


and 2 shows that the spacecraft mass capability has decreased for all launch day - trip 
time combinations. On launch days which had previously been single conic days (such 
as March 6), it is optimum to use the propellant as soon as possible after spacecraft 
separation (from the launch vehicle) for values of spacecraft mass to about 550 pounds 
mass (250 kg). It is assumed that about 1 day is required to obtain preliminary tracking 
information and proper spacecraft orientation. For the values of hyperbolic velocity 
required for this mission, the spacecraft has essentially escaped from the Earth’s 
sphere of influence after 1 day, so that the spacecraft burn takes place in heliocentric 
space. However, the spacecraft is very near to Earth, relative to heliocentric dimen- 
sions. 
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On launch days which previously had required broken plane trajectories, part of the 
spacecraft propellant is burned near the Earth, with the remainder (at least as much as 
used in the unrestricted case) used to construct essentially the same broken plane tra- 
jectory as previously. For a given launch date and trip time, the payload degradation 
resulting from the 100 pounds mass (45 kg) of spacecraft propellant is smaller when the 
unrestricted trip had been a broken plane transfer. This result is expected, since the 
additional spacecraft propellant mass is a smaller amount above the optimum for these 
cases. 

The minimum trip times for the restricted-spacecraft-propellant-mass case are 
plotted in figure 9 as a function of launch date and are illustrated for several values of 



Launch date (1972) 

Figure 9. - Minimum trip time as a function of launch date. 

Spacecraft propellant mass, 100 pounds mass (45 kg). 

spacecraft mass, as in figure 3. The increase in minimum trip time can be seen by 
comparison with figure 3. Spacecraft mass of 580 pounds mass (263 kg) or greater 
cannot be delivered when 100 pounds mass (45 kg) of spacecraft propellant must be 
carried. 

The duration of launch opportunity for the restricted case may be determined from 
figure 9 and may be compared to the results from figure 3. The number of additional 
broken plane days is decreased by at most 1 day if at all, for all values of maximum 
trip time and spacecraft mass. 

The mass of spececraft propellant utilized in a broken plane manner is shown in 
figure 10 for the launch day - trip time combinations defined in figure 9. The trajec- 
tories launched in early to middle March 1972, are single conic type for a spacecraft 
mass less than 560 pounds mass (254 kg), in that all the propellant is burned near 
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2/19 2/23 2/27 3/2 3/6 3/10 3/14 3/18 3/22 3/26 3/30 4/3 

Launch date (1972) 

Figure 10. - Mass of broken-plane-type spacecraft propellant as a function 
of launch date. Total spacecraft propellant mass (45 kg). 



Launch date (1972) 

Figure 11. - Optimum time of spacecraft correction as a func- 
tion of launch date. Spacecraft propellant mass, 100 pounds 
mass (45 kg). 
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Earth. For a spacecraft mass of 560 pounds mass (254 kg) or greater, all the tra- 
jectories are broken plane type. The mass of broken plane type spacecraft propellant 
shown in figure 10 is always at least as large as the corresponding optimum value pre- 
sented in figure 4. As in the unrestricted case, the broken plane type propellant mass 
increases to nearly 100 pounds mass (45 kg) at both ends of the launch opportunity. 

The optimum spacecraft correction times corresponding to the launch day - trip time 
combinations in figure 9 and the propellant mass in figure 10 are shown in figure 11. 
These values do not differ appreciably from those presented in figure 5. 


CONCLUDING REMARKS 

The results presented show that broken plane transfers can be used to increase the 
duration of the Pioneer- F launch opportunity. Broken plane trajectories which deliver 
the required spacecraft mass are available both before the opening and after the close 
of the single conic launch opportunity. The required trip times before the single conic 
opening are in the same range as those required for the single conic trips. However, 
the required trip times after the closing of the single conic opportunity increase rapidly 
with launch date. Therefore, the possible increase in available launch days depends 
on the maximum allowable trip time. 

To take full advantage of the broken plane trajectories, the existing spacecraft pro- 
pellant tanks must be resized to allow up to 100 additional pounds mass (45 kg) of pro- 
pellant. For best performance, the amount of spacecraft propellant loaded must be 
varied from day to day. If this is not possible, the additional 100 pounds mass (45 kg) 
is loaded and utilized on each launch day. The required trip times for this case are 
somewhat longer than if the propellant can be optimally tanked, but no more than 1 poten- 
tial launch day is lost for any value of spacecraft mass. 

The addition of spacecraft propellant will cause an increase in propellant tank mass, 
and possibly in structural mass. This mass increase has not been evaluated or con- 
sidered in this report. In addition, it may be possible in certain cases to combine a 
required midcourse correction with a planned spacecraft burn, thus reducing the total 
amount of propellant required for the two maneuvers. This may offset some, or all, of 
the additional propellant tank and structure mass. 

The results presented show that the use of broken plane trajectories gives only a 
small improvement for the Pioneer- F mission. It should be recognized, however, that 
these results are only valid for the 1972 Jupiter flyby opportunity. The advantage of 
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broken plane trajectories depends on the existing geometry, which varies from oppor- 
tunity to opportunity. However, the report serves to illustrate the tradeoffs which must 
be considered in evaluating broken plane trajectories. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 17, 1970, 

731-25. 
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APPENDIX A 


A, B, C,D 
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G* 
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SYMBOLS 


submatrices of state tran- 
sition matrix 

jet velocity, ft/sec (m/sec) 

thrust direction 

central body gravitational 
constant, ft 3 /sec 2 
(m 3 /sec 2 ) 

variational Hamiltonian, 
lbm/sec (kg/sec) 

mass, lbm (kg) 

position, ft (m) 

time, sec 

velocity, ft/sec (m/sec) 

mass flow rate, lbm/sec 
(kg/sec) 

declination, deg 

switching function 

primer vector, lbm-sec/ft 
(kg-sec/m) 

adjoint variable, lbm/ft 
(kg/m) 


$ state transition matrix 

cp right ascension, deg 

ct adjoint variable 

Subscripts: 
e earth 

f final 

H hyperbolic 

max maximum 

o initial 

2 intermediate 

Superscripts: 

derivative with respect to 
time 

— . vector 

unit vector 

~ before an impulse 

+ after an impulse 
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APPENDIX B 


MATHEMATICAL ANALYSIS 

The mathematical problem to be solved is to find optimal fixed-transfer-time im- 
pulsive transfers between two fixed positions in three-dimensional space with an 
inverse-square force field acting. This problem was studied by Pines (ref. 6), and the 
important results are repeated here for convenience. 

The equations of motion are 


^-G! ?+ C£f 

3 m 


r = V 
m = -0 


(Bla) 


(Bib) 

(Blc) 


where the controls are the thrust direction f and the flow rate 0, which are restricted 
by | f [ =1 and 0 > 0. The variational Hamiltonian is 


H = 



r" + — X* • f + jiT • V - a/3 
m 


(B2) 


and the adjoint equations are 



(B3a) 


-^91 r-s£i<r. r)r 

3 5 

r r 


a = C0r.f 


m 


To maximize H with respect to f and 0, 


(B3b) 


(B3c) 


f = X 
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and 


/3 = 0 if K = — X • f - a =— X - a < 0 
m m 

p—+°° if k > 0 

and /3 is indeterminate if #c = 0. The singular arc, k = 0, was treated by Pines but will 
not be considered here. Pines showed that impulses occur when k = 0, and furthermore, 
that k < 0 for the entire (optimal) trajectory. In addition, he showed that the value of 
ma is continuous throughout the trajectory. 

For the present problem, the impulses are supplied by the spacecraft with jet ve- 
locity C. The launch vehicle is treated by allowing variations in the initial mass and 
velocity in the transversality equation. Consider first the problem of transferring be- 
tween fixed positions ? and F^ with a fixed transfer time, and let all of the required 
velocity at F q be supplied by the launch vehicle. No impulse is required at Fj since 
the mission is a flyby with no specification on final velocity. The required AV Q is ob- 
tained by solving Lambert's problem. The transversality equation for this problem is 

^X • d V + p • dF + a dm - H dtj t - • d V + p • dF + a dm - H dt^ t - dm { = 0 

where m^ is the quantity to be maximized. The initial and final positions and transfer 
time are fixed so that 



and 


H f = H o 

The final velocity is free, so that X ^ = 0. Also, m Q (the injected spacecraft mass) is 
a function of AV = V - V g , where V 0 is the Earth's heliocentric velocity. Then dm Q 
can be expressed as 

^o = m o * d V D 


and since V Q is free, 
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■ mi i i ii in 


A = -or rn 
o o V Q o 


Also, m f is free, so that = 1; and since m 0 = m f (there is no spacecraft impulse), it 
follows that cr Q = 1. Then 



From the known values of ^ and X Q , the values of X and /i can be calculated at any 
point on the trajectory by using the state transition matrix as discussed by Pines (ref. 6) 
and presented explicitly by Danby (ref. 7). If the state transition matrix is partitioned, 

/A B 
\C D 

then X^. = AX q + B’jjL o = 0 and 11 = -B~^AX q . The values of X and "jl at any point along 
the trajectory are then obtained by application of the state transition matrix to X Q 
and £ q . 

To test the optimality of this single conic trajectory, the maximum value of X along 
the trajectory is determined. If CX mov /m < 1, then no spacecraft impulse is desired 

iii ax 

along the trajectory, and the single conic trajectory is an optimum (probably, but not 
necessarily, a global optimum). If CX m „ v /m > 1, then the single conic trajectory is 
not locally optimum. The time and position vector t 0 and r 0 , corresponding to 
are then calculated, modified slightly, and used as first estimates of the optimum time 
and position of the spacecraft impulse. Then, trajectories are constructed between F 
and Fg ancl between Fg and Fp with trip times tg - t and t^ - tg, respectively. 

The primer vector along this trajectory is 



The values of an d ^ ma y be obtained by applying the state transition matrix be- 
tween t and t 0 , and between t 0 and t.. It is shown in reference 8 that the final 
o & it 

mass is related to Fg and tg as follows: 
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(B5) 


V r 2 m f ~ ^2 " m 2 



The optimum values of r 2 and t 2 (and corresponding values of spacecraft propel- 
lant mass and direction of impulse) may be determined by using a search procedure to 
maximize m^ with respect to r” 2 and t 2 , using the gradients in equation (B 5). In this 
report, the Fletcher -Powell modified gradient method (ref. 4) was used to perform this 
optimization. 
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